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PREFACE 

The  high-frequency  electromagnetic  pulse  radiated  from  a  nuclear 
explosion  presents  a  possible  detection  mechanism  for  nuclear  explo¬ 
sions.  However,  since  receivers  for  such  signals  are  generally  at 
a  great  distance  from  the  actual  explosion,  it  is  necessary  to  con¬ 
sider  the  effects  of  ground-wave  and  sky-wave  propagation  on  the 
received  signals.  This  Memorandum  considers  these  questions  and 
helps  evaluate  the  effectiveness  of  pulse  radiation  as  a  detection 
mechanism. 

This  Memorandum  was  prepared  at  the  invitation  of  Dr.  Robert 
Frosch,  Deputy  Director  of  the  Advanced  Research  Projects  Agency, 
for  inclusion  in  the  IEEE  special  issue  on  nuclear  detection. 


V 


SUMMARY 


The  high-frequency  electromagnetic  pulse  radiation  from  nuclear 
explosions  presents  a  possible  detection  mechanism  for  low-altitude  or 
ground  bursts  as  well  as  for  bursts  in  outer  space.  This  radiation 
propagates  around  the  earth  via  ground  wave  or  ionospheric  reflection, 
and  it  is  necessary  to  determine  what  effect  this  propagation  will 
have  on  the  received  signals. 

This  Memorandum  first  considers  the  attenuation  of  the  pulse  when 
it  travels  via  ground  wave.  Effects  of  the  height  of  the  source  are 
included.  Severe  selective  attenuation  does  occur,  with  the  higher 
frequency  components  suffering  most.  The  ionospheric  waves  suffer 
relatively  little  attenuation,  but  are  strongly  dispersed,  leading  to 
limitations  on  the  pulse  length  which  can  be  resolved.  The  receiver 
bandwidth  should  be  selected  to  match  the  ionospheric  dispersion  at 
the  receiver  frequency. 

From  all  these  considerations,  it  is  clear  that  oropagation  ef¬ 
fects  will  very  strongly  modify  the  pulse  shape  and  pulse  amplitude 
relative  to  free  space. 
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I,  INTRODUCTION 

It  is  well  known  that  nuclear  explosions  produce  strong  electro¬ 
magnetic  signals. ^  It  has  been  proposed  to  use  these  signals  as  a 

(2) 

detection  mechanism  for  nuclear  bursts  in  outer  space.  The  signals 
also  constitute  a  possible  detection  mechanism  for  ground  or  low- 
altitude  bursts.  However,  receivers  for  such  signals  will  generally 
be  at  a  very  great  distance  from  the  source,  and  it  is  necessary  to 
ascertain  what  effect  electromagnetic  propagation  will  have  on  the 
received  signals. 

As  shown  in  Ref.  1,  the  signal  may  be  viewed  as  an  early  rapid 
oscillation,  and  a  later  slower  oscillation.  The  gamma  rays  from  the 
explosion  produce  Compcon  electrons,  which  ionize  the  air,  causing  the 
development  of  a  large  conducting  region.  This  charged  region  then  de¬ 
cays  slowly,  producing  electromagnetic  radiation,  with  characteristic 
wavelengths  on  the  scale  of  the  mean  dimension  of  the  ionized  region, 
provided  there  is  a  departure  from  spherical  symmetry.  In  addition, 
there  is  high-frequency  radiation  associated  with  the  initial  build-up 
phase  of  the  ionization. 

The  general  scale  of  the  ionized  region  depends  on  yield,  but  for 
sea- level  explosions  an  average  size  of  1  km  encompasses  a  very  wide 
range  of  yields.  Asymmetry  may  be  produced  by  the  proximity  of  the 
ground,  by  the  variation  of  atmospheric  density  with  altitude,  or  by 
deflection  of  the  Compton  electrons  by  the  earth's  magnetic  field. 

The  frequency  of  the  principal  mode  of  the  natural  oscillations  of  a 
sphere  of  radius  a  is  /3  c/(4rra)  —  (40/a^)  Kc,  whence  it  is  clear 
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that  the  portion  of  the  signal  arising  from  the  decay  of  the  ionization 
will  be  at  VLF  and  LF. 

The  radiation  of  such  signals  around  the  earth  has  been  considered 
rather  thoroughly  by  Wait,^^  Jfohler,^  and  their  collaborators. 

They  show  that  VLF  signals  are  so  modified  by  propagation  that  at 
sufficiently  great  distances  virtually  all  signals  look  alike. 

The  detailed  structure  of  the  source  signal  will  not  be  considered. 

However,  the  electronic  conductivity  associated  with  the  ionization 

electrons  will  decay  by  attachment  to  oxygen,  and  this  process  has  a 

characteristic  time  of  .01  M<sec  at  sea  level.  The  high-frequency  source 

will  therefore  have  a  rapid  rise,  determined  by  the  gamma- ray  production 

rate,  and  will  be  essentially  gone  in  about  1  P-sec.  Depending  on  the 

altitude  of  burst,  the  peak  electric  field  intensity  in  volts  per  meter 

3  6 

at  a  distance  D  km  from  the  burst  point  may  be  between  10  /D  and  10  /D 
for  free-space  propagation  conditions. 

High-frequency  electromagnetic  radiation  may  propagate  over  the 
earth  via  ground  wave  or  ionospheric  reflection.  Thus,  a  relatively 
close  receiver  could  observe  the  ground  wave,  and  a  more  distant  one, 
working  at  frequencies  below  the  maximum  operating  frequency  for  the 
path  between  the  burst  and  the  receiver,  could  observe  one-  or  multiple- 
hop  sky  waves. 

The  ground  wave  will  experience  severe  attenuation  as  it  propagates 
over  the  earth,  with  the  higher  frequencies  suffering  most.  The  iono¬ 
spheric  waves,  reflected  and  transmitted,  suffer  relatively  little  at¬ 
tenuation,  but  are  strongly  dispersed,  leading  to  limitations  on  the 
pulse  length  which  can  be  resolved.  Selecting  340  km  as  a  characteristic 
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distance  for  ground  waves,  and  3000  km  for  sky  waves,  the  modifications 
of  a  pulse  in  propagating  from  a  reference  distance  of  1  km  to  an 
appropriately  selected  receiver  shall  be  ascertained. 
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II.  PROPAGATION  OF  PULSES  VIA  GROUND  WAVE 

While  there  have  been  very  thorough  investigations  of  steady-state 
ground  wave  propagation,  the  work  on  pulse  propagation  is  quite  limited. 
Aside  from  the  previously  mentioned  VLF  work,  there  have  been  theoretical 
investigations  of  HF  pulse  propagation  by  Levy  and  Keller,^  and  by 
Gyunninen.  This  section  *ill  first  consider  a  steady-state  investi¬ 
gation,  and  then  will  follow  with  certain  improvements  the  method  of 
Levy  and  Keller. 

The  steady-state  theory  is  based  on  the  curves  of  attenuatijn  versus 
frequency  and  distance  presented  in  the  treatise  on  wave  propagation  by 

(7) 

Bremmer.  7  Information  has  been  extracted  from  these  curves  to  present 
Fig.  1,  which  shows  as  a  function  of  frequency  the  attenuation  below  free 
space  of  a  steady-state  signal  for  a  340  km  sea-water  path.  The  three 
curves  represent  transmitter  heights  of  0,  2,  and  4  km,  with  the  receiver 
on  the  surface. 

The  curve  for  a  transmitter  on  the  ground  has  been  read  directly 
from  Ref.  7.  The  elevated  source  curves  are  obtained  by  assuming  free- 
space  propagation  from  the  source  to  its  horizon,  then  finding  the  ground- 
wave  attenuation  for  the  remainder  of  the  path.  This  procedure  is  valid 
if  the  receiver  is  many  wavelengths  below  the  horizon,  which  it  is  for 
the  curves  as  presented. 

It  is  apparent  from  Fig.  1  that  the  attenuation  increases  rapidly 
for  frequencies  above  about  5  Me.  This  will  produce  a  high-frequency  cut¬ 
off  in  the  received  signal,  with  a  consequent  loss  in  measurement  accuracy 
for  the  rise  time  of  the  signal.  Thus,  it  the  attenuation  curves  of  Fig. 
1,  at  least  for  the  source  on  the  surface,  wern  approximated  by  a  sharp 
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Fig.  1— Attenuation  below  free  space  for  340  km  path  over  sea  water 
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cutoff  at  5  Me,  the  fastest  rise  time  which  could  be  measured  by  an  un¬ 
compensated  wideband  receiver  would  be  about  0.1  p,sec.  Any  faster  rise 
would  be  lost.  A  careful  examination  of  the  curves  of  Ref.  7  shows  that 
this  nominal  cutoff  frequency  is  inversely  proportional  to  the  cube  of 
the  distance  from  the  source  horizon  to  the  receiver.  Thus,  at  1000  km, 
the  shortest  measurable  rise  time  would  be  about  2  p,sec. 

A  representative  attenuation  for  the  significant  frequency  compo¬ 
nents  may  be  taken  as  10  db,  which  applies  at  1.7  He  for  the  ground 
source,  5.5  Me  for  the  2  km  sourca,  and  10  Me  for  the  4  km  source.  The 

amplitude  of  the  corresponding  frequency  component  at  a  distance  of  340 

-3 

km  is  then  about  10  its  amplitude  at  the  source.  If  this  is  regarded 
as  representative  of  the  peak  amplitude,  a  peak  of  1  to  1000  volts  per 
meter  may  be  observed,  with  the  lower  values  corresponding  to  low-altitude 
burs ts. 

This  steady-state  analysis  is  based  on  attenuation  versus  frequency 
curves.  If  it  is  desired  to  do  a  pulse  analysis,  it  is  necessary  in  ad¬ 
dition  to  have  phase  versus  frequency  information.  The  convolution  of 
the  source  waveform  with  the  impulse  response  of  the  ground-wave  trans¬ 
mission  system  then  yields  the  shape  of  the  received  pulse.  This  con¬ 
volution  can  be  done  numerically,  but  it  is  then  very  difficult  to  obtain 
any  generalizations.  Rather  than  this,  the  analytic  procedure  of  Levy 
and  Keller^^  will  be  followed.  Certain  improvements  have  been  developed 
in  the  procedure,  and  the  relevant  analysis  is  presented  here.  The  no¬ 
tation  has  been  slightly  modified  to  make  use  of  existing  tables  of 
functions. 

The  method  first  obtains  the  Hertz  vector  which  represents  the 
field  due  to  a  time-harmonic  soiree,  then  superposes  such  solutions  to 
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obtain  the  field  of  a  delta- function  source.  Following  Ref.  5,  suppose 
a  vertical  electric  dipole  is  located  at  the  point  r  =  p,9=0ofa 
polar  coordinate  system  with  its  origin  at  the  center  of  the  earth. 

The  fields  can  be  represented  in  terms  of  derivatives  of  a  Hertz  vector 
which  has  only  a  radial  component  r  U.  If  the  field  is  time-harmonic  with 
angula  frequency  a),  then  the  Hertz  vector  U  =  u  e  Let  e,  H>,  and 

cr  denote  dielectric  constant,  permeability,  and  conductivity,  with  a 
subscript  1  denoting  values  inside  the  earth  (r  <  a) ,  and  0=0  outside. 

The  propagation  constants  k  and  k^  are  defined  by 


r—  2  2 

k  =  iuv'e  |i  =  uj/c  ,  k^  =  a) 


(1) 


A  parameter  6  will  be  defined  by 


6  = 


2  2 
i  kJVk* 


<k  a)1/3[<k2/k2)-lf2 


(2) 


In  general,  the  conductivity  is  sufficiently  high  that  the  displacement 
current  in  the  ground  may  be  neglected.  The  parameter  6  then  simplifies 
to 


6  = 


(377  Uj.) 


1/2 


(ka) 


-5/6 


(3) 


For  sea  water,  this  is  quite  closely 

i  — 

6  =  e  ^  x  10^  (ka)  (4) 

Important  values  of  6  turn  out  to  be  near  unity.  Thus,  the  important 
values  of  ka  may  be  expected  to  be  ne ir  10^,  which  corresponds  to  a 
frequency  of  7.5  Me. 
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For  large  ka,  the  function  u(r,9)  can  be  expressed  as 


i? 


.  *  fk.l1/6  r  exp[ie{k**-T  (ka)1/3)] 

u(r>e)  ■  <%2_  l  - - 


H(r)H(p) 


(2nsin0)  a  g_j^ 


2t  -  6 

s 


(5) 


Here  the  height-gain  function  H(r)  is  equal  to 


2ni 


2tt1 

K(r)  -  Al[21/3«  3  {t8-  ^  <ka)2/3}]/Ai{2l/3a  3  T,} 


(6) 


The  Ai  or  Airy  function,  which  replaces  the  w  function  of  Levy  and 

/o\ 

Keller,  has  been  investigated  and  tabulated  for  real  arguments.  ' 
The  numbers  Tg  are  the  roots  of  the  equation 


2ni 


Ai  ' 

[V'3e  3  T  1 

_ _ s- 

Ai[21/3  e“T  T  J 


i  2 
3 


(7) 


As  6  varies,  the  first  two  roots  of  this  equation  lie  in  the  ranges 

.  T!  .  TT 

1  3  *3 

Tl  *  (.8086  to  1. 8557)e  J,  t2  =  (2.5781  to  3.2446)e  (8) 

For  the  distances  and  values  of  ka  of  interest,  the  second  and  higher 
terms  in  the  expansion  of  Eq.  (5)  have  large  negative  exponential  factors 
and  may  be  neglected  relative  to  the  first.  The  subscript  "s"  will  be 
dropped  hereafter  and  referred  to  as  t.  As  is  apparent  from  Eq.  (7), 
t  is  a  function  of  6,  and  thence,  through  Eq.  (3),  of  ka.  The  deriva¬ 
tive  of  t  with  respect  to  6  can  be  explicitly  represented  in  terns  of 


t  and  6  by 
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dr  - - 1 - 

d6  2t62  -  1 


(9) 


The  solution  of  the  pulse  problem  is  obtained  by  superposing 

•  i»lcc  t 

solutions  of  the  form  u  e  .  The  Hertz  vector  of  a  delta-function 

source  in  the  absence  of  the  earth  is  given  by 


U  »  6(t-D/c)/4nD 
o 


1 

4ttD 


c 

2tt 


dk  e 


-  ik(ct-D) 


(10) 


where  D  is  the  distance  from  source  to  receiver.  In  the  presence  of 
the  earth,  the  field  of  this  source  is 


u « r  dk 

2rr  «j 


-ikct  . 
e  u(r,0) 


(ID 


Substituting  from  Eq.  (5)  and  dropping  the  sumnation  sign,  there  results 


U  = 


c  e 


.  tt 

1  4 


2rra(2nsin6) 


1  J 

2-00 


(*  dk(ka)6  H(r)H(p)e~  ik<ct“*9>  +  i0T(ka) 


2t  -  6 


-2 


(12) 


To  this  point  the  work  of  Levy  and  Keller^  has  been  followed  di¬ 
rectly.  The  next  step  is  to  evaluate  the  integral  in  Eq.  (12)  by  the 
saddle-point  technique,  which  regards  the  integrand  as  the  product  of 
an  exponential  fa' tor  and  a  slowly  varying  factor,  finds  the  value  of 
k  for  which  the  exponent  has  a  maximum,  and  expands  to  seconu  order 
about  that  maximum.  The  resulting  integral  may  then  be  evaluated 
explicitly.  Levy  and  Keller  did  not  include  all  of  the  exponential 
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parts  of  the  integrand  in  their  evaluation,  and  consequently  their 

results  are  in  error.  Tills  error  will  be  rectified  here. 

The  height-gain  function  exhibits  exponential  behavior  for  large 

2/3  4 

values  of  its  argument.  Since  (ka)  ri.,.y  be  on  the  order  of  i. 0  ,  the 
argument  may  be  large  for  heights  of  a  few  kilometers.  The  asymptotic 
form  for  the  Airy  function  presented  in  Ref.  8  is 


_  1  3 

Ai(z)  -  z  ^  exp(-  z2) 


Thus,  the  receiver  height-gain  function  H(r)  may  be  written  as 


H(r)  -  F(r)  exp  i-j^j  j  ka  -  (ka)3}- 


with  a  corresponding  representation  for  the  source  height-gain  function 
H(p).  The  function  F(r)  is  algebraic  in  ka  for  large  values  of  ka. 

When  these  representations  are  inserted  into  the  integral  of  Kq. 
(12),  it  may  be  written  as 


2Tra(2TTsin9) 


1  J 

2  -  « 


dk(ka)6  F(r)F(p)e~  ik(ct~p/>  +  iTes  (ka> 
2t  -  6“2 


In  Eq.  (15),  the  distance  D'  and  angle  9  are  defined  as 

s 


11 


e 

s 


=  e 


PW 


riifii £i 

a 


(17) 


and  have  simple  interpretations.  The  distance  d'  is  the  shortest 

distance  from  source  to  receiver  in  free  space,  obtained  by  drawing 

straight  lines  from  the  source  to  its  horizon  and  from  the  receiver 

to  its  horizon,  then  connecting  these  by  a  circular  arc.  This  circular 

arc  has  the  gular  length  0  ,  and  represents  the  "shadow  distance"  in 

s 

which  all  tne  attenuation  below  free  space  occurs.  These  interpreta¬ 
tions  are  all  subject  to  the  approximation  that  the  source  and  re¬ 
ceiver  altitudes  are  small  compared  to  the  radius  of  the  earth,  which 
is  the  only  condition  under  which  ground-wave  propagation  can  occur. 

Let  Dg  =  a0g,  t*ie  shadow  distance. 

The  integral  of  Eq.  (15)  vanishes  for  ct  <  D',  which  is  a  necessary 
condition.  Define  the  time  T  =  t-(D'/c)  as  the  time  measured  at  the 
receiver  from  the  first  arrival  of  signal.  The  exponent  in  Eq.  (15) 
has  an  extreme  value  when  its  derivative  with  respect  to  ka  vanishes. 

The  corresponding  value  of  k  will  be  called  kQ.  The  variation  of  t 
with  ka  must  be  included,  and  leads  to  the  equation 


c'i 

D 

s 


d6 

d(ka) 


(koa) 


1/3 


(18) 


Substituting  for  the  derivatives  from  Eqs.  (3)  and  (9)  yields 


3cT 

D 

s 


(k  a) 

o 


-2/3 


M 


2t62-1J 


(19) 
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For  large  values  of  6,  the  second  term  in  the  bracket  of  Eq.  (19)  may 
be  neglected.  Since  t  then  h.;s  phase  tt/3,  kQa  must  have  phase  rr/2  to 
make  the  right  side  real.  Referring  to  Eq.  (3).  6  has  phase  n/3  when 
k^a  has  phase  tt/2.  Equation  (7)  then  indicates  that  t  and  6  continue 
to  have  phase  tt/3  as  ka  decreases,  and  for  such  values  the  bracket  on 
the  right  side  of  Eq.  (19)  also  retains  the  phase  n/3.  Therefore,  a 
new  parameter  t*  will  be  introduced,  defined  by 


T 


t 


_ UJ _ 

(2|t||6|2+1) 


(20) 


The  saddle  point  then  occurs  at 

,  T'D  3/2 

M  "  (  3CT-  )  (2 

As  ka  decreases  from  infinity  to  zero,  |t|  increases  from  0.8086  to 
1.8558.  However,  increases  more  rapidly  than  |t|,  reaching  a 
maximum  value  of  2.183  at  |t|  =  1.5,  then  decreasing  to  ?.8558. 

The  exponent  in  Eq.  (15)  is  .  expanded  to  second  order  about 
the  saddle  point  kQ.  The  factors  F(r),  F(p),  and  (ka)  '  3re  all  ap¬ 
proximated  by  their  values  at  the  saddle  point,  and  removed  from  the 
integrand.  The  remaining  integration  may  be  performed  immediately. 
After  many  algebraic  reductions,  the  expression  for  U  is 


U  = 


3/2 


-icSrcM-io, 


1/2 


F(r)F(p) 


23/2na(DDs)1/2  (2|t||6|2+  1)  t'+|  |o| 


1/2 


(22) 
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Here  F(r),  the  algebraic  height-gain  function  evaluated  at  the  saddle 


point,  is 


(\  r,,  N_3/2  .  D  t'  3/2  Xll/2  Dr'  l/2l 

PCr)  -  -MPfSij  (*-)  J 


and  there  is  a  corresponding  expression  for  F(p) . 

These  equations  are  sufficiently  complex  that  the  method  of  com¬ 
putation  requires  explanation.  First,  choose  the  distance  D  and  heights 

r-a,  p-a.  Compute  the  distances  D/  and  D  from  Eqs.  (16)  and  (17). 

s 

Select  a  sequence  of  values  of  ji|  in  the  range  .808  to  1.856,  so 
1/3  1 

chosen  that  2  j  t  j  takes  on  simple  steps.  (The  table  on  p.  477  of 
Ref.  8  gives  the  Airy  function  and  its  derivative  in  steps  of  0.1 
through  the  range,  and  this  table  is  suitable  for  interpolation.) 
Ascertain  Ai  and  Ai ‘  from  the  table,  then  compute  1 6  j  from  Eq.  (7). 
Compute  t*  from  Eq.  (20),  and  -gr^y  by  differentiation.  These  computa¬ 
tions  enable  the  factors  in  Eq.  (22)  which  only  depend  on  J t j  and  j 6  I 


to  be  obtained  once  and  for  all.  Figure  2  presents  rf  and  G( | T J )  versus 
| T | ,  where  G  is  given  by 

1/2 

G(M)  -  |6|2/{2|t|!6|2+  1}  {t'+  I  1 6 1  }  (24) 


From  1 6  j ,  |kQa|  may  be  determined  by  Eq.  (4).  Figure  3  presents  ) 6 | 
and  )kQaJ  x  10  ^versus  | t  J ,  in  which  sea-water  conductivity  has  been 


employed. 
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Fig. 2— The  functions  r'and  G  (Irl) 
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Fran  the  values  of  (k^a) ,  t7,  and  D^,  the  value  of  T  may  be  de¬ 
termined  by  Eq.  (19)  ,  thereby  establishing  the  time  after  the  arrival 
of  the  wavefront  that  corresponds  to  a  given  value  of  |t|.  This  is 
about  as  far  as  one  can  proceed  in  terms  of  universal  curves.  Tne  re¬ 
maining  factors  of  U  must  be  evaluated  separately  for  each  value  of  D, 
r,  and  p.  The  special  case  D  =  340  km,  r  =  a  has  been  worked  out  in 
detail  for  several  values  of  p,  and  the  resulting  pulse  shapes  are  pre¬ 
sented  in  Fig.  4.  The  ordinate  is  amplitude  in  arbitrary  unite,  and 
the  source  altitude  in  kilometers  is  the  labeling  parameter. 

These  curves  indicate  that  the  shortest  rise  time  observable  by 
a  wideband  receiver  at  340  km  from  a  burst  is  about  .10  p,sec  for  a 
ground  source,  .05  p,sec  for  a  burst  at  2  km  altitude,  and  .03  p,sec  for 
a  4  km  burst.  These  numbers  are  quite  consistent  with  the  interpreta¬ 
tion  of  the  steady-state  theory  given  earlier  in  this  section.  The 

/  3  .  \ 

main  factor  in  determining  the  time  of  the  peak  is  the  term  t  / 3cT J 
in  the  exponential  factor  of  Eq.  (22)  which  demonstrates  the  cubic  de¬ 
pendence  of  rise  time  on  distance.  The  drop  in  amplitude  of  the  re¬ 
ceived  pulse  as  the  source  height  increases  from  zero  to  . 75  km  occurs 
because  the  algebraic  height-gain  function  initially  decreases  with 
altitude.  The  Hertz  vector  of  the  peak  value  of  the  pulse  waveform 
decreases  as  the  inverse  fourth  power  of  the  distance  from  the  source. 

The  electric  field  associated  with  the  Hertz  vector  contains  a  factor 

2 

(kQa)  ,  and  hence  the  peak  value  decreases  as  the  inverse  tenth  power 
of  the  distance.  (See  Eqs.  (21)  and  (22).)  These  results  are  for  a 
delta- function  source.  Similar  conclusions  apply  to  a  source  waveform 


with  finite  rise  and  fall  times 


17 


9pD4!|duiv 


Fig. 4— Waveshape  of  a  delta-function  pulse  for  various  source  altitudes 
as  propagated  for  340  km  path  over  sea  water 
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It  is  apparent  from  these  considerations  that  the  selective  at¬ 
tenuation  of  high  frequencies  by  ground-wave  propagation  will  cause  a 
rapid  decrease  in  the  possibility  of  high-frequency  detection  and  wave¬ 
form  interpretation  as  the  distance  from  tie  source  increases  signifi¬ 
cantly  beyond  the  340  km  considered  here.  These  results  are  for  a  sea¬ 
water  pa?h.  Over  a  land  path,  the  attenuation  is  very  much  greater, 
and  the  possibility  of  HF  detection  very  much  less.  The  sky-wave  propa¬ 
gation  characteristics  will  now  be  investigated  to  ascertain  what  hap¬ 
pens  to  the  pulse  waveform  in  propagating  to  great  distances  via  iono¬ 


spheric  reflection. 
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III.  PROPAGATION  OF  PULSES  AROUND  THE  EARTH  VIA  THE  IONOSPHERE 


Volumes  have  been  written  on  steady-state  wave  propagation  around 
the  earth  via  ionospheric  reflection.  The  National  Bureau  of  Standards 
publishes  information  monthly  on  the  maximum  usable  frequencies  for 
various  paths.  However,  analysis  or  experiments  on  transmission  of 
the  high-frequency  tail  of  short  pulses  is  very  limited. 

Even  the  steady-state  ionospheric  theory  indicates  a  very  complex 
pattern  for  long-range  high-frequency  propagation.  Multiple  hops  and 
the  presence  or  absence  of  certain  ionospheric  layers  may  make  propaga¬ 
tion  over  long  paths  highly  erratic.  When  pulses  are  transmitted,  the 
multiple  reflections  may  cause  the  appearance  of  several  pulses  at  the 
receiver.  The  effects  to  be  considered  will  bo  analyzed  on  the  assump¬ 
tion  that  the  receiver  frequency  and  bandwidth  are  selected  so  that  the 
multiple  pulses  do  not  overlap. 

Under  conditions  such  that  D-  and  E- layer  absorption  is  negligible, 
ionospheric  propagation  at  high  frequencies  to  a  distance  of  about  3000 
km  displays  relatively  little  attenuation.  Hence,  the  amplitude  of  the 
received  signal  falls  essentially  as  inverse  distance.  The  signal  may 
typically  be  from  6  to  10  db  below  free  spa-'e.  However,  the  height  of 
reflection  is  a  function  of  frequency,  and  the  signal  components  at  dif¬ 
ferent  frequencies  experience  different  amounts  of  delay.  This  disper¬ 
sive  effect  causes  an  initially  sharp  pulse  to  be  spread  out  in  time 
and  reduced  in  amplitude.  This  change  in  pulse  shape  will  be  analyzed 
here,  under  certain  approximations  which  may  be  proved  mutually  con¬ 


sistent. 
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The  ionospheric  characteristic*  for  each  hop  will  be  approximated 
by  a  constant  attenuation  factor  K,  and  a  time  delay  T  which  is  a  func¬ 
tion  of  frequency.  The  source  pulse  will  be  represented  by  its  Fourier 
spectrum  A(f).  Then  the  signal  arriving  at  the  receiver  is  given  by 


E(t)  »  |  J  df  A(f)  ei2TTf[t-T(f)l 


(25) 


The  integral  extends  over  those  values  of  frequency  for  which  trans¬ 
mission  is  possible. 

Take  the  receiver  as  flat- topped,  of  bandwidth  B  and  center 

frequency  F.  The  bandwidth  B  will  be  taken  as  small  in  the  sense  that 

the  variation  of  the  input  spectrum  A(f)  over  the  bandwidth  B  will  be 

neglected,  and  A(f)  approximated  by  A(F).  For  a  delta- function  pulse, 

which  has  a  flat  spectrum,  this  is  no  restriction.  For  a  pulse  which 

has  a  rise  time  t^,  and  a  duration  t^ ,  the  spectrum  will  be  approximately 

constant  for  frequencies  below  will  decrease  as  1/f  for  fre- 

2 

quencies  between  l/2TTt?  and  l/2rrt^,  and  will  decrease  as  1/f  for  fre¬ 
quencies  above  l/2rrt^.  The  bandwidth  and  amplitude  variations  will 
be  shown  to  meet  the  restrictions. 

The  phase  function  fT(f)  will  be  approximated  by  its  expansion  to 
second  order  around  the  frequency  F.  Thus,  quadratic  phase  distortion 
is  included.  After  these  approximations,  the  received  signal  is  given 
by 


E(t)  = 


B 

|  A(F>.i2"F<t'T<F»J  df 


HX*-  s Fi(F>)-  i f2  3  FT<F>] 

w  oF 


B 

2 


(26) 
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The  quantity  d/dFLFT(F)]  is  the  group  time  delay  1’  ,  while  the  second 

U 

derivative  is  equal  to  the  derivative  of  group  delay  with  frequency. 
The  latter  quantity  it,  generally  called  the  dispersion. 

The  integral  in  Eq.  (26)  may  be  evaluated  in  terms  of  Fresnel  in¬ 
tegrals.  Define  the  parameters  B  and  a  and  the  variables  z  by  the 

Ox 

relations 


Bo  =  /57 | a  Tg/^F I  (27) 

a  =  B/B  (28) 

o 

z  =  B  (t-T_)  ±  a  (29) 

±  o  o 

The  interest  is  in  the  magnitude  of  E(t),  the  phase  being  relatively 
unimportant.  This  magnitude  is  given  by 

K  B  A.(F)  2  2  1/2 

— —  |_{c(v  “  c(z-); +  is(z+) "  s(z->|  J  (30) 

where  C  and  S  are  the  Fresnel  integrals  tabulated  on  pp.  321-324  of 
Ref.  8. 

The  parameter  Bq  is  a  property  of  the  transmission  medium,  while 

B  is  a  receiver  property.  Figure  5  displays  the  z-dependent  function 

of  Eq.  (30)  plotted  against  B  (t-T  )  for  <y  =  0.5,  1,  2.  As  is  clear 

o  o 

from  the  figure,  increasing  a  from  . 5  to  1  nearly  doubles  the  amplitude 

of  the  pulse  and  narrows  its  width.  However,  increasing  a  from  1  to  2 

reduces  the  amplitude  while  increasing  the  width.  The  cause  of  this 

effect  can  be  seen  from  Eq.  (26).  The  quadratic  term  in  the  phase 
2 

has  the  value  a  rr/2  at  the  edge  of  the  band.  Thus,  increasing  a  beyond 


-4-3-2-10  1  2  3  4 


B0(t-T0) 

Fig.5 — The  dispersed  pulse  for  various  bandwidths 
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1  causes  out-of-phaee  components  to  be  added  into  the  integral,  with 

the  resultant  spreading  and  amplitude  reduction.  It  may  be  deduced 

that  the  optimum  bandwidth  for  pulse  reception  is  the  valut  B  =  B^, 

and  the  effective  puisewidth  W.  which  will  be  taken  as  the  width 

between  the  hair-voltage  points,  is  approximately  1.2/Bq.  This  is 

the  narrowest  pulse  which  can  be  received. 

With  the  bandwidth  and  pulse  shape  established,  it  is  row  required 

to  evaluate  tue  delay  and  dispersion  characteristics  of  the  ionosphere. 

Very  elaborate  ray- tracing  programs  have  been  developed  to  permit  cal- 

(9) 

culation  of  phase  and  group  delay  on  a  digital  computer.  '  However, 
the  calculation  of  dispersion  from  these  programs  is  very  difficult, 
and  it  is  more  expeditious  and  sufficiently  accurate  to  vply  an  ana¬ 
lytic  procedure  developed  by  Kift.^*^ 

Kift  replaces  the  ac*-  il  ionosphere  by  an  equivalent  parabolic 
layer  and  neglects  the  earth's  magnetic  field.  The  path  of  a  ray 
through  such  a  layer  to  its  reflection  height  and  back  down  to  the 
surface  may  be  found  analytically,  and  so  may  the  travel  time.  The 
results  are  summarized  by  the  following  notation  and  equations 
6  =  launch  angle  of  elevation  of  the  ray 
hQ  =  height  of  the  lower  edge  of  the  layer 
y  =  semithickness  of  the  layer 

D  =  ground  distance  from  transmitter  to  receiver 

R  =  radius  of  earth 

r  =  R  +  h 
o  o 

f  =  critical  frequency  of  layer 


f  ■  signal  frequency 
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x  «  f/f 

o 

U  -  y/ro 

R  . 
n  ■  —  cos  6 
r 

o 

P  ■  slant  range  or  ray  path  length 

D  -  D1  +  D2  (31) 


Dj  “  2  R  M  x  tanh  1  jng)  (32) 

D2  =  2  R(cos" 1  n  -  5)  *  2  R  0  (33) 


P 


(34) 


Ionospheric  conditions  are  generally  such  that  the  height  and 
thickness  of  th'  layer  are  both  small  compared  to  the  radius  of  the 
earth.  For  long-range  propagation,  the  angle  6  is  generally  small. 
Equations  (31),  (32),  and  (33)  may  be  regarded  as  determining  the  angle 
6  at  which  a  signal  of  frequency  f  must  be  launched  if  it  is  to  return 
to  earth  at  distance  D.  There  may  be  several  such  angles,  but  only 
the  lowest  will  be  considered.  The  equations  may  be  solved  for  5, 
regarding  6,  h/ro>  and  y/*Q  as  first-order  quantities.  The  resulting 
expression  for  6  may  be  substituted  in  Eq.  (34),  which  then  giver  the 
group  time  of  flight  =  P/c  as  a  function  of  frequency.  After  many 

VJ 

simplifications,  there  results 

2r  no  tn  \  ,  /2h  x. 

TG  =  ~T  Sin  2R  +  f  i1'005  Ir)  x  Unh‘  {-£-) 


(35) 
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The  first  term  is  the  propagation  time  of  a  ray  reflected  from 
the  lower  boundary  of  the  ionosphere.  All  the  frequency  dependence 
is  in  the  second  or  penetration  term.  Equation  (35)  can  be  differ¬ 
entiated  to  obtain  the  dispersion.  There  results 


tanh 


2h  x 
o 

D 


h  x 


f 1 -4 (h  x/D)^}2 

o 


(36) 


Equations  (35)  and  (36)  permit  calculation  of  the  group  delay 
and  the  dispersion.  With  these  results,  the  optimum  bandwidth  can  be 
found  from  Eq.  (27),  and  the  shortest  observable  pulse  will  have  width 
1.2/B.  However,  Eqs.  (35)  and  (36)  involve  many  parameters,  and  a  com¬ 
plete  presentation  would  become  extremely  lengthy.  Therefore,  certain 
representative  values  will  be  chosen  for  analysis.  An  examination  of 
ionospheric  profiles  shows  that  the  choices  hQ  =  200  km.  y  =  100  km 
fits  a  wide  range  of  ionospheric  conditions.  The  normal  incidence 
critical  frequency  can  lie  between  2.5  and  15  Me  for  the  p  layer,  de¬ 
pending  on  latitude,  time  of  day,  season,  and  position  in  the  sunspot 
cycle.  Rather  than  attempt  to  represent  all  of  these,  the  sample  val¬ 
ues  5  and  10  Me  will  be  used  for  calculation.  The  distance  D  has  been 

selected  as  3000  km.  For  this  distance,  the  minimum  group  delay 
2ro  D 

-  sin  tt  =  10.22  msec,  and  the  excess  delay  AT  over  this  minimum  will 

c  <-R 

be  quoted.  The  elevation  angle  6  is  quite  small.  The  results  are 


summari  ad  in  Table  1. 
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Table  1 


DISPERSION  CHARACTERISTICS  OF  IONOSPHERIC  REFLECTION 


A.  h  =  200  km  y  =  100 

0 

km  D  =  3000 

km  f  = 

o 

5  Me 

F(Mc) 

ATg(fissc) 

bTG  2 

sT  ^sec> 

B(Kc) 

W^isee) 

6 

3.564 

1.196 

1293 

.928 

9 

8.107 

1.838 

1043 

1.15 

12 

14.64 

2.532 

889 

1.35 

15 

23.38 

3.311 

777 

1.54 

18 

34.64 

4.219 

689 

1.74 

B.  h  =  200  km  y  =  100  km  D  =  3000  km  f  =  10  Me 
o  _ _ _ o _ 


F(Mc) 

AT^Oxsec) 

bTG  2 

— -  Oxsec)2 

\ 

B(Kc) 

1  . 

W(nsec) 

12 

3.564 

.598 

1829 

.656 

15 

5.595 

.756 

1626 

.738 

18 

8.107 

.919 

1475 

.813 

1 

■ 

11.12 

1.088 

1355 

.885 

24 

J  4.64 

1.266 

1257 

.959 

27 

18.73 

1.454 

1173 

1.02 
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Even  fo  the  widest  relative  bandwidth  (1300  Kc  at  6  Me),  and  for 
an  inverse  square  law  dependence  of  spectral  amplitude  on  frequency, 
the  amplitude  variation  across  the  band  is  only  ±  23  per  cent.  This 
result  justifies  the  approximation  of  constant  spectral  amplitude. 

It  may  be  deduced  from  Table  1  that  for  a  receiver  whose  bandwiduh 
equals  the  ionospheric  dispersion  Bq  ,  the  output  pulse  will  have  a 
width  of  0.6  to  2.0  (j,sec  at  3000  km  range  when  the  source  signal  is 
a  delta  function. 

The  attenuation  of  the  sky  wave  below  free  space  at  3000  km  is 

about  6  db,  whence  the  total  amplitude  reduction  K/D  is  about  1/6000. 

-4 

From  Eq.  (30),  the  peak  signal  output  will  then  be  about  1.5  x  10 

Bq  A(F).  It  may  be  observed  from  Table  1  that  Bq  is  approximately 

-i- 

proportional  to  F  If  the  source  pulse  waveform  is  approximated  by 
an  exponential  rise  to  a  peak  value  Ep  at  0.03  p,sec,  followed  by  an 
exponential  fall  to  small  values  at  1  p,sec,  then  tl ■.«  spectrum  may  be 
approximated  by 

A(F)  -  .177  Ep/F  [l  +  (££ff  (37) 

The  spectrum  thus  has  a  transition  near  15  Me.  At  the  frequencies  of 

Table  1,  the  source  spectrum  A(F)  is  decreasing  at  a  rate  between  F  * 

-2  -3/2 

and  F  ,  so  the  overall  output  is  decreasing  at  a  rate  between  F 

•5/2  -g 

and  F  .  For  a  15  Me  receiver  the  peak  output  would  be  10  Ep  for 

a  critical  frequency  of  5  Me,  and  2  x  10  ^  Ep  for  a  critical  frequency 

3  6 

of  10  Me.  With  Ep  between  10  and  10  volts  per  meter,  the  output 
pulse  peak  value  will  be  between  1  mv/m  and  1  v/m  with  f  =  5  Me  and 
twice  that  value  with  f  *  10  Me.  For  other  values  of  frequency, 
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the  output  peek  and  pulse  width  can  be  found  from  Table  1  and  Eqs. 
(30)  and  (37). 

These  results  are  all  for  single-hop  propagation.  For  n-hop 
propagation  divide  the  value  of  D  in  Eqs.  (35)  and  (36)  by  n,  and 
then  multiply  the  resulting  time  delay  and  dispersion  by  n.  The  fre¬ 
quency  must  be  sufficiently  low  that  ionospheric  penetration  does  not 
take  place.  All  matters  involving  polarization  of  the  signals  and 
directivity  of  ths  source  have  been  omitted  from  this  analysis. 

From  all  these  considerations,  it  is  clear  that  propagation  ef¬ 
fects  will  very  strongly  modify  the  pulse  shape  and  pulse  amplitude 
relative  to  free  space.  This  paper  Indicates  the  type  and  magnitude 
of  some  of  the  effects. 
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